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I. INTRODUCTION 

The Frenkel-Kontorova (FK) modeli^ was originally 
proposed by J. Frenkel and T. Kontorova in 1938 to study 
plastic deformations and twinning 3 . It was indepen- 
dently discovered by F.C. Frank and J.H. van der Merwe 
in 1949 to study ID dislocations 4 . A rigid substrate 
is modelled by a sinusoidal potential, and film particles 
move within it making up a chain linked with harmonic 
springs of a given equilibrium length a/. The period of 
the sinusoidal potential is assumed to be a s = 1. Given 
the horizontal positions of the film particles {x{\f =1 , the 
system energy is 
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At zero temperature this model presents a phase di- 
agram in the variables Jf/J s and a/ containing many 
tongues of commensurate phases separated by infinites- 
imal gaps of incommensurate structure^Si&S. These 
pathologies are removed at finite temperature^, but a 
continuous phase transition from a solid to a liquid-like 
structure remains with T c > 0. 

The FK model has been extended both to finite- 
width stripesii and truly 2D substrates 1 ^. More real- 
istic interatomic potentials for the film interaction have 
been employed 1 ^, some of them with a phenomenological 
"temperature dependence" in order to study structural 
phase transitions 1 ^. 

In this work the FK model is modified so as the film 
layer may curve, by allowing particles to displace verti- 
cally. Most methods of analysis of the classical FK model 
are not suitable in this case. Thus, we have resorted to 
a real space renormalization group (RG) numerical tech- 
nique, inspired on the analysis of Griffiths and ChouS^ 
and the Density Matrix Renormalization Group (DMRG) 
technique^. The extension of the DMRG concepts to a 
novel area of application may be of interest to practi- 
tioners. 

Heteroepitaxial surface growth, i.e.: addition on a 
given crystalline substrate of particles of a different ma- 



terial with a given lattice mismatch, has attracted re- 
cently considerable technological interes1*iZii&. When the 
lattice parameter is larger in the film than in the sub- 
strate, 3D islands may develop, making up quantum dots. 
There are three relevant growth modes: (a) Frank-van 
der Merwe (FvdM), consisting only of flat layers; (b) 
Stranski-Krastanov (SK), in which 3D islands grow on 
top of a flat wetting layer and (c) Volmer- Weber (VW) , 
in which the 3D islands appear directly on the substrate. 

Many theoretical models have been proposed, ei- 
ther within continuous elasticity theory^ and atomistic 
models^. We will focus on the semi-atomistic FK ap- 
proach. It has been applied to study the mechanism of 
island formation both for lD2i and 2D substrates with a 
realistic interatomic film potential^. For a given num- 
ber of particles, if all monolayers cover completely the 
substrate, a FvdM growth mode is considered. If the 
particles are distributed in islands, leaving part of the 
substrate uncovered, a VW growth is assumed. When- 
ever a number of covering monolayers (i.e.: a wetting 
layer) appears with islands on top of them, the growth 
mode is SK. In these works, no vertical deformation of a 
layer is allowed: each island consists of a given number of 
completely flat layers. This assumption is questionable 
in the SK case for the reason exposed in the following 
paragraph. 

Let a s and a/ be the substrate and film equilibrium lat- 
tice parameters, with a/ > a s . Let us consider a wetting 
layer composed of flat monolayers with lattice parame- 
ters increasing upwards. The last monolayer, therefore, 
is less stressed than the first one. Thus, it is less likely 
to develop an islanding mechanism above it. The for- 
mation of a larger film surface has a certain energy cost, 
which would not have any counterpart. But this restric- 
tion disappears if all monolayers are allowed to bend. 
Consider a flat and a curved monolayer on the same sub- 
strate with the same number of atoms: the curved one 
must have larger atomic separation. Thus, the wetting 
layer height corresponds to the value where the extra en- 
ergy needed to curve the film free surface is overcome by 
the reduction of the internal elastic energy provided by 
the effective stretching of all monolayers. 

A first step in this direction is the analysis of the be- 
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haviour of a single monolayer of film which is allowed to 
curve. This single monolayer may provide a seed for the 
future formation of islands when more monolayers are 
deposited. 

The rest of this article is organized as follows. Sec- 
tion [H] discusses the vertically extended FK model. Sec- 
tion II I II provides some analytical insight into the solu- 
tions. Section llVl presents the RG-like technique we have 
employed to get the minimum energy states and sectionM 
discusses the numerical approach, along with the phase 
diagram for the model. 

II. MODEL 

The original FK model splits the energy for each con- 
figuration of the film particles into two parts: the in- 
teraction between each deposited film particle and the 
substrate is represented using a continuous sinusoidal 
substrate potential, and interaction between nearest- 
neighbor particles of the film is considered to be har- 
monic. 

This splitting is preserved in our extension, but with 
some modifications. First, film particles are allowed to 
displace in the direction orthogonal to the substrate. 
These vertical displacements should be inhibited by an 
appropriate increase in the substrate potential energy. 
Thus, the sinusoidal substrate potential must be trans- 
formed into a realistic 1+1D potential. Also, the inter- 
particle film potential must be changed in accordance 
with the substrate in order to have the same long dis- 
tance behaviour. A harmonic film potential is not realis- 
tic: when two neighbour film particles get sufficiently far 
away, the interaction should be negligible^. 

In our model the X axis shall be parallel to the sub- 
strate and the Z axis shall be orthogonal to it. Particle i 
has coordinates = (xi, Zi) and the configuration energy 
has the form: 

N N-l 

E[{n}} = E W + E KHh+i - n\) (2) 

i=l i=l 

This model shall be referred to as Vertically Extended 
Frenkel-Kontorova (VEFK) 23 . Our target in this work 
is to find the ground states, i.e.: the global minimum 
energy configurations, and their properties. 

A. Film potential 

Anharmonic film potentials for the FK model have 
been studied for a long timaH. They must fulfill the 
following conditions: (a) as we get further away from 
the minimum, it should be harder to squeeze particles 
together than to separate them and (b) asymptotic free- 
dom: particles should be free when their distance is large. 
There are many potentials fulfilling these conditions in 



the literature, and most details concerning them are ir- 
relevant for the physics of our problem. We have chosen 
the Mie potential^. 

where d = |rj+i — is the actual distance between neigh- 
boring film particles, aj is their equilibrium distance, A$ 
is a constant and fif > Vf. In this work, fj,f = 12 and 
Vj = 6, which correspond to a Lenard- Jones potential. 

We define Kf to be the spring constant of small oscilla- 
tions around the equilibrium point of the film potential. 
In terms of Kf, we find Aq — K faj/ [ijVf. 

B. Substrate potential 

The FK substrate potential is extended in the verti- 
cal (Z) direction in such a way that it reduces to its 
previous sinusoidal form when z = z eq , a given equilib- 
rium value. Away from it, the potential is modified so as 
it is energetically favourable to stay at the equilibrium 
height, but not compulsory. The equilibrium points of 
this potential must be arranged periodically at points 
(x, z) = (na s ,z eq ), where neZ and a s is the substrate 
lattice parameter, which we shall assume to be a s = 1 
throughout this article. 

The conditions on the height dependence of the sub- 
strate potential are similar to those for the film: (a) 
asymptotic freedom in the Z-axis, (b) very high barrier 
at the substrate height. We propose the following form 
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where Bo, /i s and v s are constants. The +1 term by 
the cosine function makes the potential at z = z eq zero 
on average. For z < z eq this average value is positive, 
and it is negative, although tending to zero, for z > z eq . 
This ensures that the ground state is always linked to the 
substrate. Figure 2] shows this potential for some typical 
values of the parameters. The potential may be seen to 
be large near the substrate and zero for large heights. 
There is an array of minima with their corresponding 
corridors, separated by potential barriers whose height 
decrease with z. The contour view, in the horizontal 
plane, depicts some ellipses whose centers correspond to 
these minima. We have chosen low values for the expo- 
nents: fi s — 4 and v s — 2, so as to have deeper corridors 
and a clearer phase transition. 

The substrate particles creating this potential are ar- 
ranged at coordinate points (x, z) = ((n + l/2)a s ,0), as 
is depicted in figure |2 From the figure it is obvious that 
9 S = arctan(2z e? /a s ). 
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FIG. 1: Vertically extended Frenkel-Kontorova substrate po- 
tential. 
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FIG. 2: Black and white dots represent the substrate particles 
and the substrate potential equilibrium positions respectively. 



Let us consider a film particle staying close to a mini- 
mum of the substrate potential. From the hessian matrix 
we may find the second derivative in the direction of 9 S , 
towards one of the substrate particles. If we relate this 
value to K s , a spring constant for the substrate-film in- 
teraction, we find that 



Bn = K s 



(fi s - v s ){a 2 + <jz 2 eq ) 
2fi s (n 2 + 2{i s i/ s - 2vf 



(5) 



In the rest of this paper we shall assume a s = 1 and 



z eq = 0.35. We define k = K f /K s 
energy scale we assume K s + Kf 



In order to define the 
= 1. 



III. ANALYTICAL APPROACH 

We are going to introduce an Ansatz which provides 
insight into the physics of the model. We shall consider 
all particles to have their x-coordinates in the equilibrium 
positions of the substrate potential. Furthermore, let us 
suppose that a fraction p of them leave the equilibrium 
plane z = z eq and reach a new height value z\. These 
particles shall be termed runaway particles and we shall 
assume them to be regularly spaced^. As a particular 
case, if p = 1/2 we get the zig-zag state, in which every 
other particle leaves the equilibrium plane as a way to 
relax the film springs. 
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FIG. 3: Function e(zi) (energy per particle minus the value 
for the flat configuration) within the Ansatz exposed in the 
text. The plots have different values of k, but they all take 
a/ = 1.05 and p = 1/2. (Upper) The graphs illustrate the 
behaviour of the energy around fei, where a second (rough) 
minimum appears, and kt, where the two minima become 
equal. (Lower) The graphs show the situation around fe, 
where the first (flat) minimum becomes unstable. 



Let us consider the function e(zi) to be the energy per 
particle minus the value for the flat configuration, for 
given values of p, k and aj. Figure 01 (upper) shows this 
function for five values of k, with a/ = 1.05 and p = 1/2. 
The first two cases, k < k\ and k = fei « 2.261, have 
a single minimum at z± — z eq . Therefore, the ground 
state is flat. At k > k\ a second minimum appears at 
z\ > z eq . The ground state for the third curve, with 
ki < k < kt ~ 2.905, is flat, but there is also a metastable 
rough state. At k = kt ~ 2.905 the two minima become 
equal. This point marks the rugosity transition. When 
k > kt the ground state is rough, but the flat state is still 
metastable for some range in k. Figure (lower) shows 
that this situation extends up to k = k% w 8.128, when 
the flat minimum becomes a saddle point. Beyond that 
value, for k > &2, the flat state becomes unstable, and a 
new metastable state appears, with z\ < z eq . The true 
ground state is rough and has z\ > z eq . 
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FIG. 4: Phase transition, kvsaf, for the VEFK model within 
the Ansatz exposed in section 11111 The full line labelled as 
k t denotes the phase transition. Above that line, the ground 
state is rough. Below the ki curve the ground state is flat. Be- 
tween ki and kt, although the grue ground state is flat, there 
is also a metastable rough state. Between kt and ki, the true 
ground state is rough, but the flat state is still metastable. 
Beyond ki, the ground state is rough and the flat state is 
unstable. 



Now let us consider the phase diagram, i.e.: for each 
value of ay, the value of k for which the two minima of 
the function e(zi), at z = z eq and z\ > z eq , take the same 
energy. This transition curve kt(a,f) is depicted in figure 
01 with a full line. The lower curve shows fci(a/), the 
values of k for which a second local minimum appears at 
z > z eq . Between these two curves, there is a metastable 
rough state. The upper curve shows fc2(a/), the points 
where the flat state becomes unstable. Therefore, in the 
region between the transition curve and this line, the flat 
state is metastable. 

The energy e(z\) is proportional to p. Therefore, p 
plays an irrelevant role in the determination of the tran- 
sition point. This means that, within this Ansatz, the 
phase transition is not a truly collective phenomenon. 
Since p — 1/2 is the maximum possible value of p, this 
Ansatz predicts that beyond the transition point the 
rough states will always be zig-zag states. 

The transition curve ktifif) follows with good approx- 
imation a power law kt ~ (a/ — l)' 3 with two different 
exponents. For low values of a/ we have (3 « — 1.30±0.02, 
and for high values (3 ~ —1.94 ±0.02. The exponents are 
not universal: they depend on the exponents of the sub- 
strate and film potentials. We shall compare these results 
with the numerical ones. 

Although this Ansatz is quite crude, it provides cer- 
tainly some insight into the physics of the model. The ba- 
sic framework, with a transition curve and the metasta- 
bility regions, shall remain valid in the more accurate 
approach of section [Vj We shall see that p = 1/2 is not 
the preferred state after the transition. The fraction of 
runaway particles plays a more relevant role when the 



atoms are allowed to move in the A-axis. 



IV. RENORMALIZATION GROUP 
PROCEDURE 

A. Formulation 

The main procedure involved by real space RG is to 
rewrite the system equations in terms of block variables 
without altering their form. As it has been stated, our 
RG technique is strongly inspired on the Griffiths and 
Chou 5 approach and the Density Matrix Renormaliza- 
tion Group (DMRG)i&. 

A single particle moves within our 2D substrate po- 
tential, V s , which we shall term from now on. The 
particle stays at position n, and the system energy is 
V^'(ri). Let us add a second particle on (e.g.) the right 
side, and link them with the "spring" Vf(ri,r 2 ). The 
second particle shall remain fixed at any point r 2 and we 
shall force the first particle to find its equilibrium posi- 
tion by minimizing V^{r\) + Vf{r\,r 2 ). Thus, we may 
write r\ = /i^ 2 (r 2 ) f° r a known /i<_2 function. The 
first particle will be said to be enslaved to the second one 
from now on, and function /i«_2 shall be called a slaving 
mapping. 

Now, the total energy of the system may be expressed 
as: 

V^(r 2 ) = V^(n) + V f ( ri ,r 2 ) + V s (r 2 ) 
= ^ (1) (iWr 2 )) + y/(/i^ 2 (r 2 ),r 2 ) + V s (r 2 ) (6) 

i.e.: the total energy of the two-particles block depends 
only on the position of the second particle. Thus, it may 
be considered as an effective or renormalized potential. 
The minimum value of is therefore the minimum 
possible energy for all the two-particles system. Mini- 
mizing for a single particle, we obtain the minimum for 
a two-particles chain. 

The process may, of course, be iterated until a large 
chain is represented this way. If we have V^-^On-i), 
we may add a ro-th particle and find the value of r„_i 
which minimizes I/(™ _1 ) (r ra _x) + Vf(r n -i,r n ) for each 
value of r n . This is again a slaving relation: r n _i = 
/„_i^„(r„). The (n — l)-th particle is now enslaved to 
the n-th one. Now let us define V^ n \r n ) as 

V™(r n ) = V^irn-i) + V f {r n -x,r n ) + V s (r n ) 

= V (n-1) (/»-i,-»(r„)) + V f {f n -^ n {r n ),r n ) + V s (r n ) 

(7) 

Using the same argument as before, the minimum 
value of y( n ) is the minimum possible energy for the 
whole chain. 
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B. Reconstruction Procedure 

Let us suppose that we have carried out the previous 
procedure up to a given value of n and found the series 
of potientials V"M . We may reconstruct the 

ground state in this way. The value r° which minimizes 
(r„) is the position of the n-th particle in the ground 
state. When we have r„, we may find r„_i by minimizing 

F(»- 1 )(r„_ 1 )+V f (r n _ 1) r°) (8) 

Which is equivalent to the use of the enslaving function 
r°_j = /„_!<_„ (r°). Iterating backwards we can find the 
position of the first particle: 



If the flow took place in a finite-dimensional vector 
space, this equation would be sure to have a solution 24 . 
In our problem, despite the absence of a rigorous proof, 
we have found it to exist for a wide variety of cases. 
Solving equation ^] is a hard problem unless the RG 
procedure, starting with = V s converges. This con- 
vergence has always taken place in the cases we have 
tried. 

If we are interested only in the thermodynamic limit, 
we may find the ground state from the fixed point poten- 
tial R. The procedure is identical to the one described 
in the IIV Bl section, but using always the same potential. 
If finite-size effects are thought to be important, the full 
series of potentials should be used. 



r? = A-a(/ a -3(---(/n-i-n(r° ))•••)) (9) 

and, thus, obtain the exact solution {r®, . . . , r }. 

Effectively, relation J7J) may be expanded backwards, 
substituting with its definition in terms of V^ n ^ 2 \ 

etc. Just as a check, we may evaluate V^ n '{r^) and find 

V^(r° n ) = V s (r n ) + V f (r n _ 1 ,r° n )+ 

71 — 1 11 — 2 

+ E^ ) + E t/ /^ > r °+i) ( 10 ) 

»=1 z=l 

Which is, as it was already stated, the energy of the 
full chain in the ground state. Thus, the minimization 
equation for r n within potential V^ n ' solves the problem 
for the n-particles block, which is the essence of a real 
space RG method. 



C. Renormalization Group Flow and Fixed Points 

We may consider to be flowing under RG trans- 
formations in a certain functional space. Although the 
parameter space is infinite dimensional, this poses no spe- 
cial theoretical difficulties. All potentials are periodic in 
x with period a s , since a block displacement of the full 
chain of a s does not affect the energy, but they need not 
even be continuous. We will define a fixed point in this 
RG flow through the equation: 

V^(r)=V^(r) + e (II) 

The parameter e is just the energy per particle in the 
full chain. This fixed point equation may be rewritten as 
an eigenvalue equation within a minimax algebra— for a 
target potential R{ri): 

R{n) = Vsir^ + mbxlVfiln-n-iD+Riri-x)] +e (12) 



V. NUMERICAL APPROACH 

The most straightforward strategy to obtain the 
ground state is to obtain the zero-force equations for 
each particle, solve them and choose the solution with 
lowest energy. The equation for depends both on r^_i 
and ri+i. This yields a dynamical system which, for the 
original FK model, is related to the standard map. A 
minimization algorithm must be used to choose the orbit 
which corresponds to the minimum energy. In our case, 
this strategy is not appropriate, since it is not possible 
to solve analytically the i-th zero- force equation for r^+i, 
as it is in the original FK model. Therefore, the nu- 
merical approach shall be based on global minimization 
techniques, whose problem is the existence of a complex 
landscape of local minima. 

In order to obtain a good seed configuration to insert 
in a global numerical minimization routine, we have tried 
several different algorithms: (a) Simulated annealing 25 

(b) Take a small chain and find its ground state; after 
that, add a site both at the left and the right borders and 
repeat the procedure until the desired length is reached; 

(c) A technique based on the Puncture Renormalization 
Group 26 : start with any initial configuration, pick up a 
block of sites and minimize its energy while the left and 
right tails may only move as two solid blocks; afterwards, 
the chosen block moves rightwards and leftwards, sweep- 
ing the system, and (d) The RG method described in 
section Hvl whose implementation we are going to detail. 

Method (d) has always provided the lowest energy min- 
ima, for reasons which will be discussed in the following 
paragraph. The other three methods give, in many cases, 
higher energy local minima. Method (b) is the second 
best, and its main handicap is that its solutions tend 
to have too much left-right symmetry. It is always con- 
venient to optimize the solution found using a classical 
minimization algorithm 30 . 
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A. Technical Details 

The procedure described in section llVl may not be car- 
ried out analytically. Our numerical implementation dis- 
cretizes each V^' into a regular grid of N x x N z points, 
typically 100 x 100 or 100 x 150 points. The discretiza- 
tion interval of the potential is an important issue. We 
have found the election [— a s /2, a s /2] x [z eq /2, 5z eq /2] to 
be appropriate. The X-interval is motivated by the a s 
periodicity of all potentials. In the Z-axis, the interval is 
designed to include all runaway particles. 

In their seminal work, Griffiths and Choi£ obtain the 
minimum enthalpy state instead of the minimum energy 
one. This means that, in our case, the film potential 
would be changed into 

Vf(n,r 2 ) = min V/(|ri - r 2 + m ■ e x \) (13) 

where e x is the unit vector in the A"-direction and m is 
an integer. This change is made so as the search for the 
minimum in the construction of V^ n+1 ^ from V^ 71 ' has 
a solution in the [— a s , a s ] interval of the X-axis. Other- 
wise, the search should be extended to the whole real line. 
We have found an alternative approach to this problem. 
The search is carried out on two periods in the X-axis 
on the left side of the (n + l)-th particle. Since a/ < 2a s 
in all the studied cases, this procedure is sure to find 
the solution. Our solution is always optimal within the 
subspace of coarse-grained potentials on the given grid. 

Our calculations are of two types: infinite systems, i.e.: 
systems in the thermodynamical limit, and finite-size sys- 
tems. The differences between them shall be analyzed in 
the following sections. 

For each value of a/ and k = Kt/K s the ground state 
is found. For each ground state the following observables 
are considered: (a) energy per particle e, (b) average 
height h = (z) and its standard deviation er^, known as 
roughness, and (c) effective lattice parameter a, defined 
as the average of \xi+i — sc»|, and its standard deviation 

B. Roughness transition: infinite systems 

For any given value of a/, the value k = Kf/K s = 
always results in a flat configuration, with h = z eq and 
a = a s . If we increase k, all the observables increase 
smoothly up to a certain value fc t (a/), where an abrupt 
jump takes place. 

The infinite system computation starts with the appli- 
cation of the RG flow on the substrate potential V s in 
order to find the fixed point potential R(r). Typically, 
about 10 or 20 iterations were enough to get a maximum 
error of a part in 10 6 . This fixed point potential yields a 
unique enslaving function, /. We obtain the point where 
R(r) takes its minimum value, r®, and iterate the enslav- 
ing function on it so as to reconstruct the full ground 
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FIG. 5: The phase diagram for the VEFK model in the ther- 
modynamical limit. 



state {r®, r®, ■ ■ ■ , r%} for a large value of N. If, for a 
given tolerance, there is a value No such that if n > No, 
z n « z eq , we shall assume that the ground state is flat. 
We have chosen a tolerance of 10~ 4 . 

The phase transition obtained with this criterion is 
shown in figure [S] The first part, up to w 1.075, 
may be fitted to fct(a/) ~ (a/ — with $ w —1.39, 
in reasonable agreement with the value —1.30 ± 0.02 ob- 
tained for the analytical Ansatz. For high values of af 
the transition curve is more rough, and the fit is more dif- 
ficult. Nonetheless, the exponent may be seen to increase 
in absolute value. 

The transition always takes place at a higher value of 
k than in the analytical Ansatz. This is due to the fact 
that particles were not allowed in that case to move in 
the equilibrium plane. This way to relax stress delays 
the transition in practice. 

According to the results of section ITTll for most points 
just above the transition the ground state is a zig-zag 
state, i.e.: the fraction of runaway particles is either p = 
1/2 or p = 1. 



C. Roughness transition: finite-size systems 

Finite-size systems present interesting peculiarities. 
Whenever an island develops and new film particles are 
deposited on it, the size of the island renders the thermo- 
dynamical limit irrelevant. Therefore, finite-size effects 
might play an important role in practical applictions. We 
have mapped and analyzed the roughness transition for 
a finite-size system with N = 101 particles. 

In these calculations, the intermediate potentials of the 
RG flow are stored to disk so as they may be recovered 
during the reconstruction procedure. When convergence 
is attained, the program stops calculating any more po- 
tentials and assumes the rest to be equal to the last one 
in storage. 
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FIG. 6: For 101 particles and a/ = 1.03, ground states with 
(a) k = 3.04670, (b) k = 3.04674 and (c) k = 3.04678. 

The positions of the particles always belong to the 
N x x N z grid. Therefore, their precision is low. In the 
finite-size case the positions entered a standard minimiza- 
tion algorithm. Usually, this procedure only took a very 
short time, since the configuration was nearly optimal 
and its deviation was only due to the effect of the lattice. 
Checks were periodically made to ensure that the energy 
per particle in our results were lower than in the other 
approaches described at the beginning of this section. 

Figure shows three ground states with a/ = 1.03 and 
k = 3.04670, 3.04674 and 3.04678 for TV = 101 particles. 
The first one is flat. In the second one, a single run- 
away particle appears. In the third one, two symetrically 
placed particles leave the equilibrium plane z = z eq . The 
values of er^ increase from 7 • 10~ 4 to 0.032 and 0.045 
respectively in a very small interval for k. As we increase 
k, new peaks appear steadily, in a slower way. 

Figure0shows the behaviour of the energy per particle 
e, roughness ah and deviation of the effective lattice pa- 
rameter a a for cif = 1.03 as we increase k. The transition 
is clearly seen to coincide for the three observables. 

For each value of k and a/ a certain flat domain-size 
D may be defined as the minimum distance between run- 
away particles, or between runaway particles and the ex- 
treme of the chain. Of course, a ground state is rough 
whenever D < N. 

Figure [H] shows the phase transition curve in the k — af 
plane for N — 101. This curve presents some jumps, 
e.g. at a/ = 1.0418 between k w 2.04 and 2.28, which 
are finite size effects. At this special jump, the domain 
size is D = 50, i.e.: half the system size. Thus, the first 
particle to leave the equilibrium plane is the central one. 
This behaviour is locked for a certain range of k values. 
The second jump in size is at a/ = 1.064, and there the 
domain size is D = 34, i.e.: very close to one third of the 
total lattice size. This locking is a collective phenomenon 
which does not appear in our analytical Ansatz. 

This behaviour is clearly different from the one we have 
found for infinite systems. Border effects are clearly dom- 




FIG. 7: The transition in three observables: (upper) the en- 
ergy per particle, which has a small jump; (middle) the rough- 
ness, which has a jump of 2 orders of magnitude; (bottom) the 
deviation in the effective lattice parameter, which also shows 
a jump. The roughness was the guide to find the transition. 

inant in the phase transition for finite systems. 

VI. CONCLUSIONS AND FURTHER WORK 

A modification of the classical FK model, known 
as Vertically Extended Frenkel-Kontorova (VEFK), has 
been proposed in which particles may displace vertically, 
thus allowing the film layer to bend. Physical patholo- 
gies associated with the FK model, such as the sharp 
commensurate-incommensurate transition, disappear in 
this case. This modified model is used to study a rough- 
ness transition which might be the seed for island forma- 
tion in surface growth phenomena, both with analytical 
and numerical approaches, at T = 0. The numerical 
technique is inspired on the effective potentials method 
of Chou and Griffiths^ and on real space renormaliza- 
tion group methods, specially on DMRG. The analytical 
approach shows both a transition curve and two broad 
metastability regions. 

This roughness transition presents some general fea- 
tures: (a) it starts with some runaway particles which 
leave the equilibrium plane, (b) the size of the flat re- 
gions (domains) near the transition point depends on the 
parameters af and k, and is rather sensitive to finite size 
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FIG. 8: The phase diagram for the VEFK model with N = 
101 particles. Above the line, the ground states are always 
rough, and below they are flat. The phase border may be 
seen to have jumps at certain positions, as it is explained in 
the text. 

effects. 

As discussed at the introduction, rough configurations 
may be the origin of domains which make up the char- 
acteristic isfands in SK and VW growth modes. The 
effect of many film monolayers is, of course, of great im- 
portance, and has been left for further work. A second 
monolayer may have a qualitative effect on the transition. 

All the surface growth phenomena are, of course, 
2+1D. The dimensional extension of our model involves 
some numerical difficulties, since there does not exist a 
straightforward extension of the RG approach. The per- 
formance of DMRG-like techniques decrease when the 



particle linking graph does not have a tree topology. Ex- 
tensions of the DMRG algorithm to 2D and 3D problems 
are a field of active researches and this problem might 
well pose a possible application for them. 

This VEFK model may be studied at finite tempera- 
ture with the transfer matrix formalism. The usual ar- 
guments against ID phase transitions do not apply in 
this case, since the substrate potential imposes a long 
range order. Therefore, the roughness transition is likely 
to survive up to a finite value of T. A positive result 
in this direction is the treatment of a finite temperature 
transition for FK. 

In physical terms, there should be a second flatness 
transition at even higher values of k = Kf/K s , for which 
the film gets flat again, with the effective lattice param- 
eter a = a f . In our numerical analyses this second tran- 
sition does not appear. The reason is the following: at 
very high values of k the substrate potential is almost in- 
existent, and the reference plane for the film disappears 
in practice. All configurations in which each particle is 
at a distance a/ from its neighbours, even if it is highly 
curved, have nearly the same energy. 

Some relevant features, such as a slowly changing cur- 
vature or the second flatness transition, would be met if 
a three-body angular potential was included. This ex- 
tension is physically plausible and may lead to a richer 
phase diagram. 
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